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ABSTRACT 



The force-free parameter a, also known as helicity parameter or twist pa- 
rameter, bears the same sign as the magnetic hehcity under some restrictive 
conditions. The single global value of a for a whole active region gives the degree 
of twist per unit axial length. We investigate the effect of polarimetric noise on 
the calculation of global a value and magnetic energy of an analytical bipole. 
The analytical bipole has been generated using the force-free field approximation 
with a known value of constant a and magnetic energy. The magnetic parameters 
obtained from the analytical bipole are used to generate Stokes profiles from the 
Unno-Rachkovsky solutions for polarized radiative transfer equations. Then we 
add random noise of the order of 10~^ of the continuum intensity (Ic) in these 
profiles to simulate the real profiles obtained by modern spectropolarimeters like 
Hinode (SOT/SP), SVM (USO), ASP, DLSP, POLIS, SOLIS etc. These noisy 
profiles are then inverted using a Milne-Eddington inversion code to retrieve the 
magnetic parameters. Hundred realizations of this process of adding random 
noise and polarimetric inversion is repeated to study the distribution of error in 
global OL and magnetic energy values. The results show that : (1). the sign of a is 
not influenced by polarimetric noise and very accurate values of global twist can 
be calculated, and (2). accurate estimation of magnetic energy with uncertainty 
as low as 0.5% is possible under the force-free condition. 

Subject headings: Sun: magnetic fields, Sun: photosphere 
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Introduction 



Helical structures in the solar features 



by George E. Hale in 1925 ((Hale 



1925 . 



ike sunspot whirls were reported long back 



19271 ). He found that about 80% of the sunspot 



whirls were counterclockwise in the northern hemisphere and clock wise in the southe rn 



194l|) by 



hemisphere. Later, in 1941 the result was confirmed by Richardson (iRichardsonl 
extending the investigation over four solar cycles. This hemispheric pattern was found 
to be independent of the solar cycle. Since the 90's, the subject has been rejuvenated 



and this hemispheric behaviour independent of su nspot eye 



for many of the sola. r features like active regions (ISeehafei 



Longcope et al. 



1998: 



Abramenko et al. 



1996 : iBao 
t al. 2OO3I : 



Bao fc Zhang 



e is claimed to be observed 



1990 



19981: 



Pevtsov et al 



1995 



Hagino fc Sakurai 



2OO5I), 



filaments (IMartin et a. 



Rust fc Kumai 



Rust 



1996 



1994: 



Pevtsov et a 



Pevtsov fc Longcope 



1994), coronal 



( iPevtsov et al. 



2001 



Bernasconi et al. 



2OO5I ). coronal loops 



20011). interp 



X-ray arcades (jMartin fc McAllister 



Pevtsov fc Longcope 



aneta ry magnetic clouds (IMCs) 



19961 ) and network magnetic fields 



20071) etc. 



Helicity is a physical quantity t hat measures the degree of linkages and twistedness in 



the field lines 



Berger fc Field 



19841 ). Magnetic helicity H^ is given by a volume integra l 



over the scalar product of the magnetic field B and its vector potential A (lElsasser 



19561). 



H,n= I A-BdV 

with B= Vx A. 



V 



It is well known that the vector potential A is not unique, thereby preve nting the 



calcul ation of a unique value for the magnetic helicity from the equation (1). 



Seehafer 



( 119901 ) pointed out that the helicity of magnetic field can best be characterized by the 



force-free parameter a, also known as the helicity parameter or twist parameter. The 
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force-free condition (j Chandr asekharl 1 1 9 6 ll : |Parkeiill979l ) is given as 



V X B = aB 



(2) 



Alpha is a measure of degree of twist per unit axial length (see Appendix-A for details of 
physical meaning of alpha). This is a local parameter which can vary across the field but is 
constant along the field lines. 

Researchers have claimed to have deter mined the sign of ma gnetic helicity on the 



photosphere by calculat 



ing alpha, e.g 



< az > = < Jz/Bz > (jPevtsov et al. 



flPevtsov et al. 



19951 ). averaged alpha e.g. 



19941 ) with current density = (V x B)^. Som e 



authors have used current he licity density He = Bz ■ Jz and aav (IBao fc Zhang 



Hagino fc Sakurai 



Burnett e et al 



1998 



2004, 



(120041 ) and 



2005r). A good corr elation was found between af,est and («^) by 



LekaetaL 



(119961 ). But the sign of magnetic helicity cannot be 



inferred from the force-free parameter a under all conditions (See Appendix B). 

It is well known that the reliable measurements of vector magnetic fields are needed to 
study various important parameters like electric currents in the active regions, magnetic 
energy dissipation during flares, field geometry of sunspots, magnetic twist etc. The 



study of error p ropagation from po. 



very important (ILites &: Skumanich 



arizat i on measurements to y e ctor field parameters is 



1985 



Klimchuk et al 



1992) 



Klimchuk et al. 



1l992h 



have studied the effects of realistic errors e.g., due to random polarization noise, crosstalk 
between different polarization signals, systematic polarization bias and seeing induced 
crosstalk etc. on known magnetic fields. They derived analytical expressions for how these 
errors produce errors in the estimation of magnetic energy (calculated from virial theorem). 
However, they simulated these effects for magnetographs which sample polarization at few 
fixed wavelength positions in line wings. It is well known that such observations lead to 
systematic under-estim ation of field strength and also suffer from magneto-optical effects 



(IWest fc Hagyard 



19831 ). Whereas in our analysis, we simulate the effect of polarimetric 
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noise on field parameters as deduced by full Stokes inversion. The details are discussed in 
the section 6. 



Pevtsov et al.l (119951 ) found large variations in the global a values from repeated 
observations of the same active regions. It is important to model the measurement 
uncertainties before looking for physical explanations for such a scatter. 



In a study by 



Hagyard fc Pevtsovl (Il999l ) the noise levels in the observed fields were 



analyzed, but a quantitative relationship between the uncertainties in fields and the 
uncertainties in global a value were not established. They could only determine the extent 
to which the incremental introduction of noise affects the observed value of alpha. However, 
for the proper tracking of error propagation, we need to start with ideal data devoid of 
noise and with known values of a and magnetic energy. We follow the latter approach in 
our present analysis. 

Here, we estimate the accuracy in the calculation of the a parameter and the magnetic 
energy due to different noise levels in the spectro-polarimetric profiles. Modern instruments 
measure the full Stokes polarization parameters within the line profile. Basically there 
are two types of sp ectro-polarime t ers : (i) Spectrograph based e.g.. Advanced Stokes 



Polarimeter (ASP : 



Keller et al. 



111992^ 



Elmore et al. 



Povell fll995h : 



19921) 1 , Zuri c h Imaging Polari r neter 



Stenflol fll996h : 



Stenflo fc Kellei] fll997h). THEM 



(ZIMPOL : 



(jArnaud et a. 



Keller et al. 



1119981)), SOLIS - Vector Spectro-Magnetograph fV SM : Ijones et al 



(120031 )). Polarimetric Littrow Spectr ograph ('POLIS : 



Diffraction Limited Spectro-polarimeter 



Hinode (SOT/SP : 



Tsuneta et al. 



(DLSP : 



Schmidt et al 



Sankarasubramanian et al 



(12004 



S-MTR 



1120021); 



12003(1), 



20061)), 



(120081 ) )■ etc. and (ii) F ilter-based e.g., Irn aging Vector 



Magnetograph (IVM) at Mees Solar Observatory, Hawaii (IMickev et al. 



Magnetograph at Udaipur Solar Observatory (SVM-USO) (lGosain et al 



1996), Solar Vector 



2004, 



20061 ) etc. 



Earlier magnetographs like Crimea (IStepanov &: Severny 



1962), MSEC ( iHagvard et al. 
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vector magnetograph (IStaude et al 



1982|), HSP flMickevl ll985h. QAQ flMakita et al.lllQSSh. HSOS fAi fc Hul ll986h. Potsdam 



19911), SFT fISakurai et al 



19951 ) etc. were mostly 



based on polarization measurements at a few wavelength positions in the line wings and 
hence su bject to Zeeman satura t ion effects as well as magneto-optical effects like Faraday 



rotation (IWest fc Hagyard 



1983 



Hagyard et al 



2000) . 



The magnetic field vector deduced from Stokes profiles by modern techniques are alra ost 



free 



200l|). 



rom such effects ( jSkumanich Sz Lites 



1987 



S'anchez Almeida 



1998 



Socas-Navarro 



This paper serves three purposes. First, we estimate the error in the calculation of 
field strength, inclination and azimuth and thereafter in the calculation of the vector field 
components B^, By and B^. Second, we estimate the error in the determination of global 
a due to noise in polarimetric profiles constructed from the analytical vector field data. 
Third, we estimate the error in the calculation of magnetic energy derived using virial 
theorem, due to polarimetric noise. 

In the next section (section 2) we discuss a direct method for calculation of a single 
global a for an active region. In section 3, we describe the method of simulating an 
analytical bipole field. Section 4 contains the analysis and the results. Error estimation in 
global a is given in section 5. In section 6 we discuss the process of estimating the error in 
the virial magnetic energy. Section 7 deals with discussion and conclusions. 



2. A direct method for calculation of global a 



Taking the z-component of magnetic field, from the force-free field equation (2) a can 
be written as, 

(V X B), 

a = (3) 



B. 
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For a least squares minimization, we should have 



[a — a.. 



minimum 



or. 



an 



a 



(4) 



where a is the local value at each pixel, ag is the global value of a for the complete 
active region and N is total number of pixels. 

Since eqn.(4) will lead to singularities at the neutral lines where approaches 0, therefore 
the next moment of minimization, 



- ag)'B' 



minimum 



should be used. From eqn. (5) we have 

_d_ 

dag 

which leads to the following result. 



an 



- agYBl) = 



^ ax ay ' ^ 



(5) 



(6) 



(7) 



This formula gives a sin gle global value of a in a sunspot and is the same as aav of 



Hagino fc Sakurail (120041 ). We prefer t his direct way of ob taining global a which is 



different from the method discussed in 



Pevtsov et al. 



(119951 ) for determining a^est- The 



main differences are : (1). the singularities at neutral line are automatically avoided in 
our method by using the second moment of minimization and (2). the computation of 



Hagino fc Sakurai 



consta nt a force-free fields for different test values of a is not required. 
(2004) used a different parameter a^av to avoid the effect of Faraday rotation in sunspot 



umbrae. However, modern inversion techniques using complete Stokes profiles are free of 
this problem. 

It must be noted that one can generate different values of ag using higher moments of 
minimization, e.g., by weighting with B^, with n=3, 5, 7, ... etc. The higher moments 
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will be more sensitive to spatial variation of B^. Suc h large and complex variation of B 



is found generally in flare productive active regions (lAmbastha et al. 



1996 



Hagyard et al. 



1993 



Wang et al. 



19991 ). Thus we can try to use higher order ag as a global index for 



predicting the flare productivity in active regions. 



Finally, to compute ag we need all the three components of magnetic field which 
is obtained from the measurements of vector magnetograms. However, here we use the 
analytically generated bipole, as discussed in the following section, with known values of all 
the magnetic parameters to investigate the effect of polarimetric noise. 



3. Generation of theoretical bipole 



We use the analytic, non-potential force- free fields of the form derived by iLowl (119821 ). 
These fields describe an isolated bipolar magnetic region which is obtained by introducing 
currents into a potential field structure. This potential field is produced by an infinite 
straight line current running along the intersection of the planes y = and z = -a, where 
negative sign denotes planes below the photosphere z = 0. At the photosphere (z = 0), the 
field has the following form : 

B, = ^cos0(r) (8) 

= ^ cos0(r) - sin0(r) (10) 

where Bq is the magnitude of the field at origin and r"^ = x"^ + y"^ + . The function 
0(r) is a free generating function related to the force-free parameter a (see eqn (2)) by 

a = -— (11) 
ar 
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which determines the current structure and hence the amount and location of shear present 
in the region. By choosing (f){r) = constant = 7r/2 we can obtain a simple potential 
(current-free, a = 0) field produced by the infinite line current lying outside the domain. 
Steeper gradient of 0(r) results in a more sheared (non-potential) field. 

In equation (11) the sign on the right hand side is taken positive in the paper by Low 
(1982) which is a typing mistake (confirmed by B. C. Low, private com naunication) . We 



menti on this here to avoid carrying forward of this typo as was done in 



Wilkinson et al. 



mm . 



A grid of 100 x 100 pixels was selected for calculating the field components. The 
magnitude of field strength at origin has been taken as lOOOG and the value of 'a' is taken 
as 15 pixels (below the photosphere, z = 0). 

The simulated field components with corresponding contours are shown in the figure 1. 



Wilkinson et al. 



(119891 )) for the generation of 



Here we use the following function (e.g., 
the field components (B^, B^^, B^) : 

0(r) = r<3a (12) 

= I , r>3a (13) 

Results for the fields generated by different 0(r) are quantitatively similar. In this way we 
generate a set of vector fields with known values of a. 

Most of the time one of the bipoles of a sunspot observed on the Sun is compact 
(leading) and the other (following) is comparatively diffuse. Observations of compact pole 
gives half of the total flux of the sunspot and is mostly used for analysis. One can derive 
the twist present in the sunspot using one compact pole of the bipolar sunspot for constant 
a. Thus we have selected a single polarity of the analytical bipole as shown in figure [1] to 
calculate the twist. 
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Fine structure in real sunspots is difficult to model. Our analysis applies to the large 
scale patterns of the magnetic field regardless of fine structure. 

All the following sections discuss the analysis and results obtained. 



4. Profile generation from the analytical data and inversion 



Using the analytical bipole method (ILow 



19821 ) the non-potential force-free field 



components B^;, Bj^ & B^ in a plane have been generated and are given as in equations 
(8), (9) & (10). We have shown B^., Bj^, & B^ maps (generated on a grid of 100 x 100 
pixels) in figured] From these components we have derived magnetic field strength (B), 
inclination (7) and azimuth (^ : free from 180° ambiguity). In order to simulate the effect 
of typical polarimetric noise in actual solar observations on magnetic field measurements 
and study the error in the calculation of a and magnetic energy, we have generated the 



synthetic Stokes profiles for each B, 7 and ^ i n a grid of 



He-Line Information Extractor "HELIX" code (iLagg et al. 



00 X 100 pixels, using the 



2004j ). This code is a Stokes 



inversion code based on fitting the observed Stokes profiles wi th synthetic ones obtained 



by Unno-Rachkovsky solutions (lUnno 



1956 



Rachkowsky 



19671 ) to the polarized radiative 



transfer equations (RTE) under t 



( iLandolfi fc Landi Degl'Innocenti 



le as sumption of Milne-Eddington (ME) atmosphere 



19821 ) and local t her mo dynamical equilibrium (LTE). 



However, one can also use this code for generating synthetic Stokes profiles for an input 
model atmosphere. The synthetic profiles are functions of magnetic field strength (B), 
inclination (7), azimuth (^), line of sight velocity (fLos), Doppler width {voopp), damping 
constant (F), ratio of the center to continuum opacity [rjo), slope of the source function 
(Sgrad) and the source function (Sq) at r = 0. The filling factor is taken as unity. In 
our profile synthesis only magnetic field parameters B, 7, C, are varied while other model 
parameters are kept same for all pixels. The typical values of other thermodynamical 
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parameters are given in table [H We use the same parameters for all pixels. Further, 
all the physical parameters at each pixel are taken to be constant in the line forming 
region. Howe ver, one must remember that real sola r observations have often Stokes V area 



asymmetries (jSolanki 



1989 



Khomenko et al. 



20051) result of vertical magnetic and 



velocity field gradients present in the line forming region. This has not been taken into 
account in our simulations. 

A set of Stokes profiles with 0.5% and 2.0% noise for a pixel is shown in figure [2l 

The wavelength grid used for generating synthetic spectral profiles is same as that 
of Hinode (SOT/SP) data which are as follows : start wavelength of 6300.89 A, spectral 
sampling 21.5 mA/pixel, and 112 spectral samples. We add normally distributed random 
noise of different levels in the synthetic Stokes profiles. Typical noise levels in Stokes profiles 



obtained by Hinode fSQT/S 



intensity, Ic fllchimoto et al. 



) no rmal mode scan are of the order of 10 of the continuum 



20081 ). We add random noise of 0.5 % of the continuum 



intensity Ic to the polarimetric profiles. In addition, we also study the effect of adding a 
noise of 2.0% level to Stokes profiles as a worst case scenario. We add 100 realizations of 
the noise of the orders mentioned above to each pixel and invert the corresponding 100 
noisy profiles using the "HELIX" code. 

The guess parameters to initialize the inversion are generated by perturbing known 
values of B, 7 and ^ by 10%. Thus after inverting 100 times we get 100 sets of B, 7 & ,^ 
maps for the input B, 7 & ^ values from bipole data. In this way we estimate the spread 
in the derived field values for various field strengths, inclinations etc. First, the inversion 
is done without adding any noise in the profiles to check the accuracy of inversion process. 
We get the results retrieved in this process which are very similar to that of the initial 
analytical ones. The scatter plot of input field strength, inclination, azimuth against the 
corresponding retrieved strength, inclination, azimuth after noise addition and inversion is 
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shown in figure [3] (upper panel). Typical B^, Bj^ & B^ maps with different noise levels are 
shown in the lower panel. As the noise increases B^^, B^, & B^ maps become more grainy. 

From the plots shown in figure [3] we can see that the error in t he field strength for a giveri 



noise level decreases for strong fields. This is similar to results of 



Venkatakrishnan fc Gary 



(119891 ). As the noise increases in the profiles, error in deriving the field strength increases. 
We find that the error in the field strength determination is ~ 15% for 0.5% noise and 
~ 25% for 2% noise in the profiles. Inclination shows more noise near 0° & 180° than at 
~ 90°. The error is less even for large noisy profiles for the inclination angles between 
~ 50° — 130°. The reason for this may be understood in the following way. Linear 
polarization is weaker near and 180° inclinations and is therefore more affected by the 
noise. The azimuth determination has inherent 180° ambiguity due to insensitivity of 
Zeeman effect to orientation of transverse fields. Thus in order to compare the input and 
output azimuths we resolve this ambiguity in ^^ut by comparing it with i.e., the value of 
^out which makes acute angle with C,in has been taken as correct. We can see azimuth values 
after resolving the ambiguity in this way show good correlation with input azimuth values. 
Some scatter is due to the points where ambiguity was not resolved due to 90° difference in 

First, the ag was calculated from the vector field components derived from the noise 
free profiles to verify the method of calculating global alpha and also the inversion process. 
We have used the single polarity to calculate global alpha present in sunspot as discussed in 
section 3. We retrieved the same value of ag as calculated using the initial analytical field 
components. From the figure |H we can see that the effect of noise on the field components 
is not much for the case of 0.5% noise but as the noise in the profiles is increased to 2.0%, 
the field components specially transverse fields show more uncertainty. The vertical field is 
comparatively less affected with noise. The scatter plot in figure |H shows that the inversion 
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gives good correlation to the actual field values. The points with large scatter are due to 
poor "signal to noise" ratio in the simulated profiles. The mean percentage error in the 
further discussions is given in terms of weighted average of error. 

5. Estimation of the error in the calculation of global alpha (ag) 

We calculate the percentage error in global alpha each time after getting the inverted 
results, for both the cases when 0.5% and 2.0% (of Ic) noise is added in the profiles, by the 
following relation : 



The histogram of the results obtained is shown in figure O 

First, we inverted the profiles without adding any noise and calculated ag from 
retrieved results to compare it with the 'true' ag calculated from the analytically generated 
vector field components. We get less than 0.002% difference in the both ag values. 

For the case of 0.5% noise in polarimetric profiles we get a mean error of 0.3% in the 
calculation of ag and error is never more than 1%. Thus the calculation of ag is almost free 
from the effect of noise in this case. Hence, by using data from Hinode (SOT/SP), one can 
derive the accurate value of twist present in a sunspot. 

If 2.0% noise is present in the polarization, then maximum ~ 5% error is obtained. 
Weighted average shows only 1% error. Thus the estimation of alpha is not influenced very 
much even from the data obtained with old and ground based magnet ographs. In any event 
it is unlikely that a realistic error will be large enough to create a change in the sign of ag. 




(14) 



where a* is calculated global alpha and ag is the analytical global alpha. 
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6. Estimation of the error in the calculation of magnetic energy (E^) 



The magn etic energy has 



virial theorem (j Chandr asekhar 



jeen c alculated using virial theorem. One form of the general 



196ll ) states that for a force-free magnetic field, the magnetic 



energy contained in a volume V is given by a surface integral over the boundary surface S, 

(B-r)B]-n dS (15) 



-Ls^dV = — [ l-Bh 
8tt An / ^2 



where r is the position vector relative to an arbitrary origin, and h is the normal vector 
at surface. Let us adopt Cartesian coordinates, taking as z=0 plane for photosphere. This 
assumption is reasonable because the size of sunspots are very small compared to the radius 
of the Sun. If we make the further reasonable assumption that the magnetic field strength 
decreases with distance more rapidly than r ~'^/^ whereas a poi nt dipole field falls off as r~^. 



197J) 



then the equation (15) can be simplified to (IMolodensky 

1 If 

—B'^dxdydz = ^ {xB^ + yBy)Bzdxdy 



(16) 



where x and y are the horizontal spatial coordinates. B^;, & B^ are the vector magnetic 
field components. This equation (16) is referred as the "magnetic virial theorem". 



Thus magnetic energy of an active region can be calculated simply by 



substituting the 



1982 



1985 



deriv ed vector field components into the surface integral of equation (16) (ILow 
19891 ). Magnetic field should be solenoidal and force-free as is the case for our analytical 
field. So the energy integral is independent of choice of the origin. 

If all the above conditions are satisfied then the remaining source of uncertainty in the 
magnetic energy estimation is the errors in the vector field measurements themselves. So, 
before the virial theorem can be meaningfully applied to the Sun, it is necessary first to 
understand how the errors in the vector field measurements produce errors in the calculated 
magnetic energies. 
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Earlier, the efforts were made to estimate the errors (iGary et aLlll987l : iKhmchuk et al. 



19921) for magnet ographs hke Marshall Space Flight Center (MSFC) magnetograph. 



Gary et al. 



3) 



(119871 ) constructed a potential field from MSFC data and computed its virial 



magnetic energy. Then, they modified the vector field components by introducing random 
error s in B^, By and and recomputed the energy. They found the two energies differ by 
11%. iKlimchuk et al.l (jl992[ ) approached the problem differently. They introduced errors in 
the polarization measurements from which the field is derived instead of introducing errors 
to magnetic fields directly. In this way they were able to approximate reality, more closely 
and were able to include certain type of err ors such as crosstalk which were beyond the 



scope of the treatment by 



Gary et al. 



(119871 ). They found that the energy uncertainties are 
likely to exceed 20% for the observations made with the vector magnetographs present at 
that time (e.g. MSFC). 



Klimchuk et al 



(I1992I ) except that we 



Here, our approach is very similar to that of 
consider full Stokes profile measurements to derive the magnetic fields like in the most of 
the recent vector magnetographs e.g., Hinode (SOT/SP), SVM-USO etc. as mentioned 
earlier. We begin with an analytical field, determine polarization signal as explained in 
earlier parts, introduce the random noise of certain known levels (0.5% & 2.0% of Ic) in the 
polarization profiles, infer an 'observed' magnetic field after doing the inversion of the noisy 
profiles, compute an 'observed' magnetic energy from the 'observed' field and then compare 
this energy with the energy of the 'true' magnetic field. The percentage error is calculated 
from the following expression: 



AEr. 



■(%) 



E*-E„ 



X 100 



17) 



E E 

where is 'observed' energy and is 'true' energy. All the above processes have been 
described in detail in section 4. 



Figure [6] shows the uncertainty estimated in the calculation of magnetic energy in two 
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cases when error in the polarimetric profiles is 0.5% and 2.0% of Ic. Needless to say, we first 
checked the procedure by calculating the magnetic energy from the vector fields derived 
from inverted results with no noise in the profiles. We found the same energy as calculated 
from the initial analytical fields. 

We can see that the magnetic energy can be calculated with a very good accuracy when 
less noise is present in the polarization as is observed in the modern telescopes like Hinode 
(SOT/SP) for which very small (of the order of 10~^ of Ic) noise is expected in profiles. We 
find that a mean of 0.5% and maximum up to 2% error is possible in the calculation of 
magnetic energy with such data. So, the magnetic energy calculated from the Hinode data 
will be very accurate provided the force-free field condition is satisfied. 

The error in the determination of magnetic energy increases for larger levels of noise. 
In the case of high noise in profiles (e.g. 2.0% of Ic) the energy estimation is very much 
vulnerable to the inaccuracies of the field values. We replaced the inverted value of the field 
parameters with the analytical value wherever the inverted values deviated by more than 
50% of the 'true' values. We then get the result shown in the right panel. We can see that 
the error is very small even in this case. The mean value of error is ~ 0.7%. 

7. Discussion and Conclusions 

We have discussed the direct method of estimating ag from vector magnetograms 
using the 2"*^ moment of minimization. The higher order moments also hold promise for 
generating an index for predicting the flare productivity in active regions. 

The global value of twist of an active region can be measured with a very good accuracy 
by calculating ag. Accurate value of twist can be obtained even if one polarity of a bipole 
is observed. 
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The magnetic energy calculation is very accurate as seen from our results. Very less 
error (approximately 0.5%) is seen in magnetic energy with 0.5% noise in the profiles. Thus 
we conclude that the magnetic energy can be estimated with very good accuracy using 
the data obtained from modern telescopes like Hinode (SOT/SP). This gives us the means 
to look for magnetic energy changes released in weak C-class flares which release radiant 
energy of the order of 10^° ergs (see Appendix-C), thereby improving the statistics. 

These energy estimates are however subject to the condition that the photospheric 
magnetic field is force-free, a condition which is not always met with. We must then obtain 
the energy estimates using vector ma gnetograms observe d at higher atmospheric layers 



where the magnetic field is force-free (IMetcalf et al. 



19951). 



The 180° azimuthal ambiguity (AA) is another source of error for determining 
parameters like ag and magnetic free energy in real sunspot observations. The smaller the 
polarimetric noise, the smaller is the uncertainty in azimuth determination, thereby allowing 
us to extend the range of the acute angle method used in our analysis. On the other hand it 
is difficult to predict the level of uncertainty produced by AA. Influence of AA is felt more 
at highly sheared regions which will anyway deviate from the global alpha value. Thus, 
avoiding such pixels will improve determination of Ug. Magnetic energy calculation at such 
pixels could be done by comparing energy estimates obtained by 'flipping' the azimuths 
and choosing the mean of the smallest and the largest estimate of the energy. Here we 
assume that half the number of pixels have the true azimuth. This is the best one can do 
for a problem that really has no theoretical solution allowed by the Zeeman effect (but see 



also, 



Metcalf et al. 



chromospheric chirality (ILopez Ariste et al. 



(12006) and references therein) . Observational techniques such as use of 



2006 



Martin et al. 



2008 



Tiwari et al. 



2008 



or 



use of magnetograms observed from different viewing angles could perhaps resolve the AA . 



Patches of both signs of alpha can be present in a single sunspot (jPevtsov et al. 



1994 : 
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Hagino &: Sakurail l2004j ). In those cases the physical meaning of ag becomes unclear. 
Efforts are needed to understand the origin of such complex variation of a in a sunspot. 
Real sunspots show filamentary structures. If this structure is accompanied by local 
variations of a, then does the global a result from correlations in the local a values? Or, 
are the small scale variations due to a turbulent cascade from the large scale features? The 
answers to these questions are beyond the scope of our present study. Modeling sunspots 
with such complex fine structures is a great challenge. However, we plan to address the 
question of fine structure of twists in real sunspots observed from HINODE (SOT/SP), in 
our forthcoming study. 

For the present, we demonstrate that the global twist present in an active region can 
be accurately measured without ambiguity in its sign. Furthermore, the high accuracy of 
magnetic energy estimation that can be obtained using data from modern instruments will 
improve the probability for detecting the flare related changes in the magnetic energy of 
active regions. 
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A. Physical meaning of force- free pcirameter a 

(Derived from the discussions with Professor Eugene N. Parker during his visit to 
Udaipur Solar Observatory) 

Taking surface integral on both sides of eqn. (2), we get 

a f dS-B = J dS-V xB 

— ^ dl -B {from Stokes theorem) (Al) 



or, 



or, 



« - (A2) 



In the cylindrical coordinate we can write eqn.(A2) as 



a 



2B. 



wBz 

where z and w are axial and radial distances from origin, respectively 

The equation of field hues in cylindrical coordinates is given as : 

Bz _ B(f, 
dz zudcf) 

B^ wdcf) 
Bz dz 

Using eqns. (A3) & (A5), we get 



(A3) 



(A4) 



(A5) 



a = 2^ (A6) 
dz 

From equation (A6) it is clear that the a gives twice the degree of twist per unit 
axial length. If we take one complete rotation of flux tube i.e., = 27r, and loop length 
A ~ lO^meters, then 
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2 X 27r , 
a = (A7) 



comes out of the order of approximately 10 ^ per meter. 



B. Correlation between sign of magnetic helicity and that of a 

Eqn. (2) can be written as 

V X B = a{V X A) 

= V X (aA) (Bl) 

giving vector potential in terms of scalar potential as 

A = Ba-i + (B2) 

which is valid only for constant a. 

Using this relation in eqn.(l), we get magnetic helicity as 

Hm = J {Ba-^ + V0) • B dV 

= j B'^a-^dV + j dV (B3) 

Second term in the right hand side of eqn. (B3) can be written as, 

J{B-V)(f)dV ^ jv- ((/)B) dV 

= J{(j)B)-ndS (B4) 

(from Gauss Divergence Theorem) which is equal to zero for a closed volume where magnetic 
field does not cross the volume boundary (n • B = 0) provided that 4> remains finite on the 
surface. Therefore, we get magnetic helicity in terms of a as 

Hm^ J B^a-^ dV (B5) 
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which shows that the force free parameter a has the same sign as that of the magnetic 
hehcity. However, if n ■ B 7^ 0, then the contribution of the second term in eqn. (B3) 
remains unspecified. Thus it is not correct to use alpha to determine the sign of magnetic 
hehcity for the half space above the photosphere since n ■ B 7^ at the photosphere. 



C. Estimate of energy release in different classes of X-ray flares: 



With the simpl ifying 



duration of 16 min (Drake 



assum ption that all classes of soft X-ray flares have a typical 



197ll ). we can see that the energy released in the different classes 



of flares will be proportional to their pe ak power. Since X- class flares typically release 



radiant energy of the order of 10'^^ ergs (lEmslie et al 



20051 ). therefore M-class, C-class, 
B-class and A-class flares will release radiant energy of the order of respectively 10^^, 10^°, 
10^9 and 1028 ergs. 



Table 1. Model parameters for generating synthetic profile 



Model Parameter 


Value 


Doppler velocity, vi^os {ms~^) 





Doppler width, voopp (mA) 


20 


Ratio of center to continuum opacity, r]o 


20 


Source function. So 


0.001 


Slope of the source function, Sgj-^d 


1.0 


Damping constant, F 


1.4 
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Fig. 1. — Contours of the field components overlaid on their gray-scale images. The contour 
levels are 100, 500 and 800 G of magnetic fields. The red and blue contours denote the 
'positive' and 'negative' polarities, respectively. The green box in shows the area which 
is selected for the calculation of global a. For details see the text. 
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Fig. 2. — Example of Stokes profiles with 0.5% (left column) and 2.0% (right column) noise 
along with fitted profiles. The input parameters for the associated pixel are as follows : field 
strength^ 861 G, inclination^lOl", azimuth=19°. The corresponding output parameters are 
850 G, 101°, 19° for 0.5% noise and 874 G, 99°, 19° for 2.0% noise. 
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Fig. 3. — Scatter plot (upper panel) between the field strength, inclination and azimuth 
before and after inversion with 0.5% (1** column) and 2.0% {2"-'^ column) noises in the profiles. 
The lower panel shows the images of vector fields B^;, By & before (1** row) and after 
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Fig. 4. — Scatter plot between Bj., By & before and after inversion without noise (I''* 
row) and with adding noise in the profiles: 2"-'^ row with 0.5% noise and 3^'^ row with 2.0% 
noise (of Ic) in the polarimetric profiles. 
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Fig. 5. — Histogram of the percentage error in calculation of ag with 0.5% and 2.0% noise 

(of Ic) in polarimetric profiles, respectively. 
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Fig. 6. — Histogram of the percentage error in calculation of magnetic energy when 0.5% 
and 2.0% noise (of Ic) is present in polarimetric profiles, respectively. 
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